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A theory of the magnetic field driven (semi-) metal- 
insulator phase transition is developed for planar systems 
with a low density of carriers and a linear (i.e., relativistic 
like) dispersion relation for low energy quasiparticles. The 
general structure of the phase diagram of the theory with re- 
spect to the coupling constant, the chemical potential and 
temperature is derived in two cases, with and without an 
external magnetic field. The conductivity and resistivity as 
functions of temperature and magnetic field are studied in de- 
tail. An exact relation for the value of the "offset" magnetic 
field B c , determining the threshold for the realization of the 
phase transition at zero temperature, is established. The the- 
ory is applied to the description of a recently observed phase 
transition induced by a magnetic field in highly oriented py- 
rolytic graphite. 

71.30.+h, 



I. INTRODUCTION 

Although during recent years there has been important 
progress in understanding non-Fermi liquid dynamics in 
dimensions D > 1, an understanding of them is still very 
far from being complete. It is rather clear that non- 
Fermi liquid behavior yields examples of sophisticated 
nonperturbative dynamics which should be described by 
advanced methods of quantum field theory. 

It was recognized rather long ago that relativistic field 
models can serve as effective theories for the descrip- 
tion of long wavelength excitations in condensed matter 
systems (for a review, see Ref. ^). In particular, they 
can be applied to a wide class of (quasi-)planar systems. 
In this case, the corresponding relativistic theories are 
(2 + l)-dimensional, i.e., they are formulated in (2 + 1)- 
dimensional Minkowski space with two space like and 
one time like coordinates. It is important that amongst 
these condensed matter systems are such as high-T c su- 
perconductors and carbon-based materials (for a list of 
papers using relativistic field approach to these systems 



see Refs. 

In this paper, we will develop a consistent approach 
to studying these systems by making use Qi-iiSQ called 
reduced (3 + l)-dimensional gauge theories These 
theories will share the following common feature. Their 
gauge fields (e.g., the electromagnetic field) responsible 
for interparticle interaction would be able to propagate 
in a three-dimensional bulk, while fermion fields (e.g., de- 
scribing electron- and hole-type quasiparticles) would be 
localized on two-dimensional planes. A typical example 
of a condensed matter system of this type is graphite. It 
has been known for a long time that fermionie-nuasipar- 
ticles in graphite are nearly two dimensional.113 In addi- 
tion, graphite is a semimetal whose low-energy quasipar- 
ticles have nearly linearpdjapersion law (just like massless 
relativistic particles). E3£3'Q The Coulomb interaction be- 
tween quasiparticles is provided by gauge fields which, 
unlike the quasiparticles themselves, are three dimen- 
sional in nature. 

Recently, the dynamics of reduced QED was studied in 
Refs. |l^|Pl| . In those papers, purely relativistic theories 
were considered: in particular the velocities of both mass- 
less fermions and photons were equal to the speed of light 
c. In realistic condensed matter systems, the Fermi ve- 
locity of gapless fermions Vp is of course much less than c. 
This in turn implies that the static Coulomb forces pro- 
vide the dominant interactions of fermions. This feature 
makes quite a difference in the analysis. 

In this paper we will describe such "realistic" reduced 
gauge theories with and without an external magnetic 
field perpendicular to the basal plane. We are particu- 
larly interested in the possibility of a spontaneous gen- 
eration of a gap in the one-quasiparticle spectrum. This 
might be viewed as a (semi)metal-insulator phase tran- 
sition. The influence of the magnetic field, as would be- 
come clear in a moment, is very powerful in driving (or 
"catalyzing" ) this type of transitions. 

The phenomenon of the magnetic catalysis of dynam- 
ical symmetry breaking was established as a universal 
phenomenon in a wide class of (2 + 1)- and (3 + 1)- 
dimensional relativistic models in Refs. [l4 15 (for ear- 



1 



the mechanism of the magnetic catalysis. Moreover, the 
exact relation for B c will be pointed out. It is: 



\eBri = 



Nt 



(1) 



where Nf is the number of fermion species ( "flavors" ) and 
n is a charge density of carriers (Nf = 2 in graphite). 
While the existence of this exact relation is noticeable 
in itself, its experimental verification would be a critical 
check of the validity of the magnetic catalysis scenario in 
highly oriented pyrolytic graphite. 

The paper is organized as follows. In Sec. |l| general 
features of the model (reduced QED) are described. In 
Sec. Ill we analyze the gap equation and establish the 
phase diag ram in reduced QED without magnetic field. 
In Sec. IV the gap equation in reduced QED with an ex- 
ternal magnetic field is studied. The resistivity and con- 
ductivity in this system are studied in detail in Sec. [v|. 
Sec. VI is devoted to the interpretation of the experimen- 



lier consideration of dynamical symmetry breaking in a 
magnetic field see Refs. |l6| , |l7| ). 

The general result states that a constant magnetic field 
leads to the generation of a fermion dynamical mass (a 
gap in a one-particle energy spectrum) even at the weak- 
est attractive interaction between fermions. The essence 
of this effect is the dimensional reduction D — ► D — 2 
in the dynamics of fermion pairing in a magnetic field. 
At weak coupling, this dynamics is dominated by the 
lowest Landpirlevel (LLL) which is essentially (D — 2)- 
dimensionalJlJ'Ej The applications of this effect-have been 
considered both in condensed matter physicscffl and cos- 
mology (for reviews see Ref. |l8| ). 

The main motivation of the present study was the ex- 
perimental data reported in Refs. [l9|-|2ll and their in- 
terpretation (based on the phenomenon of the magnetic 
catalysis) suggested in Ref. |[ It was observed in those 
experiments that samples of highly oriented pyrolytic 
graphite in an external magnetic field show a qualita- 
tive change of their resistivity as a function of temper- 
ature, that was interpreted as a metal-insulator phase 
transition. The effect is clearly seen only for a magnetic 
field perpendicular to the basal plane, suggesting that 
the orbital motion of quasiparticles is responsible for the 
change of the conductivity dependence. 

The suggestion of Ref. || was that this phenomenon 
can be a manifestation of the magnetic catalysis, when a 
dynamical gap, connected with a quasiparticle-hole pair- 
ing, is generated in a magnetic field. In this paper, we 
will develop a detailed theory of the magnetic-field-driven 
metal-insulator phase transition in planar systems, based 
on reduced QED. The general structure of the phase dia- 
gram of such systems will be described in two cases, with 
and without an external magnetic field. The behavior 
of the electric conductivity (resistivity) in these systems 
will be described in detail. This will allow us to conclude 
that, in the presence of a magnetic field, the generation 
of a dynamical gap in planar systems can indeed manifest 
itself as a metal-insulator phase transition in the behav- 
ior of the resistivity p(T, B) as a function of the magnetic 
field and temperature. 

It will be also shown that there exist clearly distin- 
guishable signatures of different types of the phase tran- 
sition. While the resistivity p(T) is a smooth function at 
the critical point T — T c in the case of a non-mean-field 
second-order phase transition, there are a discontinuity 
and a kink in p(T) at T = T c in the cases of the first-order 
and mean-field phase transitions, respectively. The con- 
clusion of the present analysis concerning the possibility 
of the realization of the scenario of the magnetic catalysis 
in highly oriented pyrolytic graphite is quite positive. 

One of the central results of this paper is an explana- 
tion of the exisjteu.ce of an "offset" field B c observed in the 
experimentsJiatil As we will discuss in detail in Sec. VI, 
the value B c determines the threshold for the generation 
of a dynamical gap at zero temperature: it happens only 
if B > B c . It is remarkable that, as will be shown in FIG. 1. The schematic lattice structure of a single layer of 

Sec. IV, the existence of B c is a robust consequence of graphite. 



tal data in highly oriented pyrolytic graphite. In Sec. [VTj, 
we summarize the results of this work. There are also 
three Appendices. The symmetry of (2 + l)-dimensional 
fermions is considered in Appendix |A|. A derivation of 
the polarization function and the gap equation in reduced 
QED is done in Appendix |b]. In Appendix ^|, the effec- 
tive potential for reduced QED with a nonzero chemical 
potential is derived. 



II. MODEL 

In this section, we describe the general features of 
the model. As mentioned in Sec. ||, the main assump- 
tion of the reduced dynamics of the planar systems is 
that the fermionic quasiparticles are confined to a plane, 
while the gauge fields are free to propagate in the three- 
dimensional bulk. 

A similar setting was recently studied in a class of rela- 
tivistic models in Refs. Here, however, we consider 

a strongly nonrelativistic model (with the Fermi velocity 
vp being much less than the speed of light) which could 
be applied to realistic planar condensed matter systems 
such as highly oriented pyrolytic graphite; see Fig. |l]. 
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The spatial coordinates on the plane (e.g., a single 
layer of graphite) are denoted by r = (x, y). The orthog- 
onal direction is labeled by the z coordinate. Thus, the 
most general bulk spatial vector is given by R = [x, y, z). 

The lagrangian density of the electromagnetic field (in 
the bulk) is given by 

C em = ^- (e E 2 -— B 2 ) -A oP +-A-j, (2) 
8tt \ Po J c 

where £o is the dielectric constant, po is the magnetic 
permeability, Aq and A are the scalar and vector poten- 
tials. The electric and magnetic fields are 



E = -VA - -d t A, 
c 

B = V x A. 



(3) 
(4) 



The interacting terms, with the quasiparticle charge den- 
sity p and current j, were also included in the lagrangian 
density in Eq. (|^) . Now, the lagrangian density of quasi- 
particles themselves (defined only on the plane) reads 



C = v F V(t,r) 



h°{dt + i[J>) . x„ . 2 q \ lT ,, ( - 
- 17 d x - ij d y ^{t,r), 



vf 



(5) 



where ^(t, r) is a four-component spinor, ^> = ^>^°, and 
the 4x4 Dirac 7-matrices furnish a reducible represen- 
tation of the Cliffo*dr-(pirac) algebra in 2 + 1 dimensions 
(see Appendix |A]) In order to describe the situation 
with a finite "residual" density of carriers, here the chem- 
ical potential p, connected with the electric charge, was 
introduced. 

We will consider the case when the fermion fields carry 
an additional, "flavor" , index i = 1,2,..., Nt (in the ex- 
ample of graphite, Nf = 2L see Refs. Then, the sym- 
metry of the lagrangian (||) is U{2Nf) (see Appendix |a|). 

In the case of minimal coupling of the electromagnetic 
field, the quasiparticle charge density and current take 
the following explicit forms: 



The bare potential Uo(t, \r\) takes the following simple 
form: 

„ , M , e 2 5(t) f d 2 k .2tt e 2 S(t) . . 

U (t, \r\) = — ^ / -— ^ exp(ifc • rW = — ^. (11) 



so J (2tt) 5 



Note, however, that in many cases of interest (e.g., in 
the case of a finite temperature and/or a finite density 
and/or a nonzero magnetic field), the polarization effects 
may considerably modify this bare Coulomb potential. 
Thus, the interaction should rather be given by 

U(jti M ) = f! / ^ / gg«Kpj-tot + iM (12) 



£0 J 2tt J 2tt |fc| +n(w,|fc|) 

where the polarization function II (w, |fc|) is proportional 
(with a factor of 27r/£o) to the time component of the 
photon polarization tensor. 

Adding a mass (gap) term Aotptp into the action 
( |io| ) would reduce the U(2Nf) symmetry down to the 
U(N f ) x U(Nf) (see Appendix |A|). Therefore the dy- 
namical generation of a fermion gap (connected with a 
quasiparticle-hole pairing) will lead to the spontaneous 
breakdown of the U(2N f ) down to the U(N f ) x U(N f )M 
Our goal is the description of the flavor phase transition 
connected with generating the gap. We will consider the 
dynamics both with and without an external magnetic 
field. 



III. GAP EQUATION. ZERO MAGNETIC FIELD 

In this section we will describe the dynamics of the 
generation of a gap connected with a quasiparticle-hole 
pairing provided by the interaction (|12|) in the case of 
the zero external magnetic field. We will begin by calcu- 
lating the polarization function n(w,|fc|). Actually, we 
will calculate (and use in the gap equation) 11(0, \k\), 
i.e., the polarization function in instantaneous approx- 
imation. The r eliabil ity of this approximation will be 
discussed in Sec. Ill E . 



j x {t,R) = ei^^rh^r)^), 
j y (t,R) =ev F ^(t,r)^(t,r)5(z), 

i«M) = o. 



(6) 
(7) 
(8) 
(9) 



Proceeding as in Ref. [lO, the initial action can be reduced 
to the plane. Then, neglecting relativistic corrections of 
order (vf/c) 2 , we are left with the following action of 
interacting quasiparticles: 



Son — I dtd 2 r£ (t,r) 



dt I dt' I d 2 r I d 2 f 



x *(t,r)~f°y{t,r)U (t-t',\r-f<\) 
x ^(t',r')7 *(t',r'). 



(10) 



A. Polarization function 

The one-loop polarization function at finite tempera- 
ture and finite chemical potential is given by the following 
integral representation (see Appendix 0): 



11(0, fc) = 



2Te N 



£oVp 



dx 



In 2 cosh 



Rx+ M 
T 



2TR r 



tanh 



R x + n 

2T 



+ (fX -> -fl) 



(13) 



where R x — y v 2 F k 2 x(l — x) + (p), At(p) is the 
fermion gap, and T is temperature. Notice that the gap 
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is a dynamic al qua ntity, determined from a gap equa- 
tion (see Sec. [II B below), and therefore it can depend 
on both temperature and chemical potential. Note that 
throughout this paper we work in the vacuum in which 
the fermion gap is positive. 

At /i = (zero density) and T = 0, the polarization 
function becomes 



11(0, fc) 



e 2 N f 
£qV 2 f 



A 



v 2 F k 2 



2v F \k\ 



4A§ v F \k\ 
arctan ■ 



2An 



(14) 



At nonzero density and T = 0, the function in Eq. (|Tj 
reduces to 



11(0, fc) = 



2e 2 N } 



— for |fc| < fc* 

£oV F 



(15) 



n( ,fc) = ^| M | 



£ov F 



1 - 



k 2 - k 2 

m 



v 2 F k 2 



4Aq(m) 



4//uf|AH, 



x arctan ■ 



v F \ k 2 - k 2 t 



2/i 



for \k\ > fc* 



(16) 



where fc* = 2^ fi 2 — Ag(/x) jv F is proportional to the 
square root of quasiparticle density at T = 0, see Eq. ( ]77| ) 
below. As is easy to check, this polarization function has 
a very strong dependence on momentum. Indeed, while 
11(0, k) remains constant for small momenta, |fc| < A;*, 
its value drops considerably for |fc| > fc*. In the case 
of a small density of carriers, i.e., n ~ k 2 -C (A /v F ) 2 , 
this momentum dependence is particularly strong. As 
is clear from Eq. (15), for small momenta, the po- 
larization function 11(0, k) is equal to the Debye mass 
Mb and could be quite large. At the same time, 
the function 11(0, k) at intermediate values of the mo- 
menta, | A; | ~ y/k*A()(n)/vF, is smaller than Md by 
about a factor of A (n) / \J [i 2 — A\ (//) , i.e., 11(0, fc) ~ 

Md^H 2 - A 2 {h)/A (h). Finally, for |fc| > A Q (/i)/w F , 
the polarization tensor approaches the following asymp- 
tote: 



n(0,fc) 



ne 2 Nf 

4:£QV F 



\k\. 



(17) 



This observation is quite important for the proper anal- 
ysis of the pairing dynamics between electron and hole 
types of quasiparticles leading to a possible dynamical 
generation of a gap. As we shall see below, it is in fact 
the region of momenta \k\ 3> Ao(p)/v F that dominates 
in such a dynamics. This in particular implies that the 
one-loop approximation with free gapless fermions (when 
both the gap and the wave function renormalization are 
neglected) is a reliable approximation for the polariza- 
tion function in the gap equation, at least for large Nf. 
It could work reasonably well evenrfor smaller values of 
Nf of order 1, say, 2 as in graphite.Ea 



B. Dynamical gap at T = and /j, = 

Now, let us study a possibility of spontaneous genera- 
tion of a dynamical gap in the one-particle spectrum of 
quasiparticles. We begin by considering the gap equation 
for the zero-density case and zero-temperature case. Its 
explicit form reads [see Eq. (B8) in Appendix^ 



A p = A 



qdqAgK, (p, q) 
vV + (A q /v F f 



(18) 



where the approximate expression for the kernel K,(p,q) 
is given by 



and 



K (p, q) 



A 



6(p-q) , e(q-p) 



2(e v F 



»2V>/4) 



(19) 



(20) 



It is well known,Elil0 that the approximation with the 
one-loop polarization function in the kernel of the gap 
equation (the so called improved rainbow approximation) 
is reliable for large values of the number of fermion flavors 
Nf. Here, however, we will also consider the values of 
Nf of order 1, say, 2 as in graphite. It is reasonable 
to assume that this approximation works qualitatively 
(although apparently not always quantitatively) even for 
these values of Nf, providing a general insight into the 
nonperturbative dynamics of spontaneous generation of 
a gap. The analysis done-jn (2 + l)-dimensional QED 
supports this assumptionE3 

In the most important region of momenta \ k\ 3> Aq/v f 
where the pairing dynamics dominates (see below), the 
only role of the term (A q /v F ) 2 in the denominator of the 
integrand on the right hand side of Eq. ( |l8| ) is to provide 
a cutoff in the infrared region. Therefore one can drop 
this term, introducing instead the explicit infrared cutoff 
Aq/v f in the integral. This is the essence of the so called 
bifurcation approximation. As a result, we arrive at the 
following equation: 




^A n 



Aq/df P 



^A 

q " 



(21) 



Here we also introduced a finite ultraviolet cutoff A. In 
a condensed matter system, it could be taken of order 
it /a where a is a characteristic lattice size (for example, 

o 

a = 2A6A for graphite). An alternative, equally good, 
estimate of A is related to the size of the energy band 
A = t/v F where t = 2.4 eV in the example of graphite. 

The last integral equation is equivalent to the differen- 
tial equation, 



p 2 A" + 2pAi + AA P = 0, 



(22) 



with the boundary conditions: 



4 



p 2 A' I 



PJ\p=A 



= 0, 

= 0. 



(23) 
(24) 



The solution compatible with the infrared boundary con- 
dition (E3|) reads 



A„ 



A, 



3/2 



sin((5) A 



; sin — 

pv F \ 2 



1 pv F 
— In — — 
A 



, (25) 



where S = arctan V4A — 1. Notice that Ao satisfies the 
relation Ao = A p= a /i>f- The ultraviolet boundary con- 
dition (^4|) imposes another restriction, 



V4A- 1 v F A 
in — h 2o = 7r. 

2 An 



(26) 



As is clear from this equation, a meaningful solution for 
the dynamical gap, satisfying the constraint Ao <SC Avp , 
exists only for A > 1/4. In the nearcritical region, i.e., 
when V4A — 1 is small, the gap reads 



An. ~ Av p exp 



2tt 



V4A - 1 



(27) 



The condition A > 1/4 gives the critical line in the plane 
(g, Nf), where the dimensionless coupling constant is g = 
e 2 /e v F : 



9c 



8 - nNf ' 



(28) 



which means that, in absence of an external magnetic 
field, a dynamical gap is generated only if g > g cr . 
In the example of graphite, the number of "flavors" is 
equal 2. Thus, the estimate of its critical coupling gives 
g cr « 2.33. We emphasize that it is just an estimate 
obtained in the leading order in 1/Nf in instantaneous 
approximation. For Nf — 2 as in graphite, both 1/Nf 
corrections and improving the instantaneous approxima- 
tion can ce rtainl y vary this value (see a discussion in the 
end of Sec. IeTFe]). 



If highly oriented pyrolytic graphite is a semimetal in 
absence of an external magnetic field, it is clear that its 
effective coupling g e ff (defined, for example, at the en- 
ergy scale below which the Dirac type effective action 
provides an appropriate description of the quasiparticle 
dynamics) is smaller than g cr . Indeed, if the interaction 
were stronger than this, the ground state rearrangement 
(from a semimetal to an insulator state), caused by the 
particle-hope pairing, could not be prevented. 

Let us now discuss the self-consistency of our assump- 
tion that the region of momenta |fc| 3> Ao/vf is mostly 
responsible for the generation of a "small" gap A <C Av F 
in the nearcritical limit. The point is that in this regime 
the logarithm ln(At\p/Ao) ~ 27r/\/4A — 1 is large. On 
the other hand, the behavior of the integrand on the right 
hand side of Eq. ( |l8| ) is smooth as q — > 0. The smooth 
behavior of the integrand in the infrared region implies 



that the region < q < Aq/v f is too small to generate 
the large logarithm ln(Ai)j?/Ao). This logarithm [and 
therefore the essential singularity in expression (27)] is 
generated in the large region Aq/vf C 5 « A. A vari- 
ation of the kernel in the infrared region can at most 
change the overall coefficient in the expression for the 
gap. 

At this point, we would like to mention that the di- 
mensionless coupling constant in the problem at hand is 
g = e 2 /€qV F . In the gap equation, g has to be considered 
as the bare coupling constant and its value can be large. 
As was shown in Rcf. ^, in the absence of a dynamical 
gap, the corresponding renormalization group (running) 
coupling runs logarithmically to a trivial fixed point in 
infrared. In the presence of the gap, such running should 
stop at the energy scale of order Ao . This means that the 
nonperturbative dynamics shifts the zero infrared fixed 
point to a finite value. 



C. Dynamical gap at T / or ^ / 

Up to now, we have considered the case with the zero 
density and zero temperature. It is clear that the criti- 
cal value of the coupling constant should be larger than 
g cr ~ 2.33 if a nonzero density (and/or finite tempera- 
ture) are taken into account. Indeed, with increasing the 
charge density of carriers or the value of the temperature, 
the screening effects get stronger and the quasiparticle in- 
teractions get weaker. In addition, the pairing between 
quasiparticles in the two adjacent bands separated by the 
dynamical gap gradually becomes less efficient. The lat- 
ter could be clearly seen by comparing the energy gain 
from creating a gap in the spectrum and the energy loss 
of pushing up the energy of all the states in the band 
above the gap. Both effects work against the formation 
of a gap. Thus, after reaching some critical value, the fi- 
nite density or temperature effects will be so strong that 
dynamical generation of a gap will be impossible. 

When the chemical potential is smaller than the gap, 
the dynamics of the zero temperature model remains un- 
changed. Thus for all values of fi < Ao = Ao(/i)| A ,=o J the 
exact solution for the dynamical gap is the same. In our 
approximation, it is given by Eq. (|27|). In order to con- 
sider the possibility of a nontrivial solution satisfying the 
condition /1 > Ao(/i), we consider an approximate gap 
equation following from Eq. ( |l7j ) and Eqs. (B8), ( p39| ) 
taken in the limit T 



0: 




(29) 



where the infrared cutoff e is given by a larger value of 
Aq(h)/vf or sj — A\{[i)/vf- By making use of the 
same method as before, we straightforwardly derive two 
branches of the solution: 



Ao(ji) ~ Ao, 



(30) 



5 



for [i < V2Aq (here we took into account that Ao(/x) 
A for fi < A ), and 



A (aO 



A 2 



(31) 




valid for A < |^| < V^Ao- While the first branch of 
the solution in Eq. ( |30| ) describes a gap that is essen- 
tially unchanged with n, the value of the gap along the 
second branch of the solution in Eq. ( |3~T| ) increases with 
the chemical potential. For the values of the chemical 
potential in the range Ao < /i < \/2Ao, both branches 
of the solution coexist. The first branch corresponds to 
a locally stable solution (i.e., to a local minimum of the 
effective potential), while the other one — to an unstable 
solution (i.e., to a local maximum of the effective poten- 
tial). In addition, there is always a trivial solution which 
corresponds to an extremum of the effective potential at 
the origin. When both nontrivial solutions (|3C]) and ( |3~l| ) 
coexist, the extremum at the origin should be a mini- 
mum. This follows from a simple consideration of the 
topology of the effective potential. 



Sigma 

- . _>s. /// 

-0.04 
-0.08 
-0.12 
-0.16 

FIG. 2. The effective potential of the composite field a for 
a few different values of the chemical potential: /1 = (solid 
line), jj, = Ao (dash-dotted line), [i = fi c ~ 1.19Ao (dashed 
line), fi — 1.3Ao (dash-dot-dotted line), and /i = \/2Ao (dot- 
ted line). In calculation, we used Ao = 1. The values of the 
potential are given in units of NfA^/vp, and values of the 
composite field a — in units of Nf(Ao/vp) 3 ^ 2 \/ r \. 

Of course, the analysis of the gap equation alone would 
not be sufficient to prove most of the above statements. 
To support them, we derived the effective potential V(a) 
as a function of the composite field a = —{tjjijj} in Ap- 
pendix This potential is graphically shown in Fig. ^ 
for a few different values of the chemical potential. As is 
clear from the figure, we have a typical realization of the 
first order phase transition. 

Our analysis of the effective potential also allows to 
determine the critical value of the chemical potential: 

A 



(2- V2)V3 

Av F 
(2- V2)V3 



exp 



2tt 



V4A- 1 



(32) 



When the chemical potential increases from \i = fi c — 
to fi = fi c + 0, the value of the gap drops from A ~ A 
down to A = 0. 

Similarly, at fj, = 0, we could derive the value of the 
critical temperature. It also appears to be of the same 
order as Aq: 



Ao 
2 



exp 



2tt 



V4A- 1 



(33) 



Unlike the case with the chemical potential, the phase 
transition in temperature is of the second order. This 
follows both from the existence of a single-branch solu- 
tion to the gap equation and from a direct study of the 
effective potential. 



D. Reduced QED3+1 vs. conventional QED 



2+1 



Before concluding this section, it is instructive to 
compare the reduced dynamics with the "conventional" 
QED 2 +i . The gap equations in these two models are sim- 
ilar, but the interaction potentials, arc .slightly different. 



Instead of expression (12), one ha: 



C/ 3 d(i,rO 



ft 
to 



2^ 



d 2 k exp(— iujt + ik ■ r) 
(2ir) 2 fc2 + n( w Jfc|) ' 



(34) 



where is the (dimensionful) coupling constant in 
QED2+1 and, in the relevant region of momenta A/vp <C 
|fc| <§; Nfe^/eoVp, the polarization function here is es- 
sentially the same as in Eq. (|l3|), except that the di- 
mensionless coupling constant (t/£qVf is replaced by the 



dimensionful eye^vp. Comparing expressions ( |34| ) and 
(|l2|), one can see that the only difference between them 
is in the appearance of the term k 2 , instead |fc|, in the 
denominator of the former. This point makes quite a 
difference. On the one hand, it provides a dynamical ul- 
traviolet cutoff ~ Nfe 2 /eoVp in the gap equation and, 
on the other hand, since this term is suppressed in the 
region |fc| <C N /e 2 ./ 'sqVp, it is irrelevant for generating 
the gap. This implies reducing screening of Coulomb like 
interactions in QED2+1 as compared to the reduced dy- 
namics. Let us consider this point in more detail. It is 
easy to find that the dynamical gap in QED2+1 is 



ye 3 

£0 



exp 



2tt 



V4A1 



(35) 



where A3 = 2/nNf. Since this solution exists when A3 > 
1/4, it implies that the critical value of Nf is equal to 
8/tt 2.55. The same critical value for Nf was obtained 
in Ref. g 

Now, notice that the parameter A3 coincides with A in 
Eq. ( p0|) only in the limit e 2 — * 00. Thus, the reduced 
dynamics becomes equivalent to QED2+1 dynamics only 
in the maximally strong coupling limit, with e 2 — * 00. 
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Therefore, we conclude that there are important similar- 
ities and important differences between the dynamics in 
QED2+1 and reduced QED3+1. Both dynamics are inti- 
mately connected with long range Coulomb like interac- 
tions. On the other hand, since QED2+1 is superrenor- 
malizable (and therefore asymptotically free) theory, its 
nonperturbative interactions are dynamically cut off at 
the scale ~ Nfe 2 /eoVp in ultraviolet. Also, its dynamics 
is more efficient in generating a dynamical gap. Indeed, 
it corresponds to the dynamics in the reduced QED3 + i 
when the coupling constant e of the latter goes to infinity. 



This feature has 
reduced QED3 + i 



een already established in relativistic 
3 with vp = c. 



E. Beyond instantaneous approximation 

Let us now turn to the discussion of the reliability of 
the instantaneous approximation for the gap equation. 
In this approximation, the frequency dependence in the 
photon propagator is neglected. While it is certainly jus- 
tified for its free (kinetic) part, it is not immediately clear 
how good it is for the polarization function. 

Keeping the frequency dependence, the gap equation 
at zero T and fi takes th e fo llowing form in E uclidean 
space [compare with Eq. (B7) in Appendix B2]: 



dcjcPk 



a s e f dud 
£0 J (2tt 



A(w,fc) 



) 2 u> 2 + v 2 F \k\ 2 + A{uj,ky 
1 



\p- k\ + U(Cl - uJ,p- k) 
where the approximate polarization function is 



2 N f \kf 



4e Ju; 2 + v 2 F \k\ 2 



(36) 



(37) 



In the instantaneous approximation, we used lo 2 — > 
in the polarization function. Thus, the strength of the 
interaction was somewhat underestimated in such an ap- 
proximation. This in turn implies that g cr and 1/N^ r 
are smaller than their values obtained in the instan- 
taneous approximation, i.e., g cr < 4/(8 — irNf) and 
Nf > 8/tt ss 2.55 in the leading order in 1/N f . 

One should remember however that 1/Nf corrections 
could be relevant for the values of Nf of order 1. Us- 
ing the argument in Ref. 25, one may expect the vari- 
ations up to 50 percent in the value of Nj r . Therefore 
the value N^ r = 2.55, obtained here in leading order in 
1/Nf, should be considered just as a useful estimate. 



IV. GAP EQUATION. NONZERO MAGNETIC 
FIELD 

The main goal of this paper is a description of the 
magnetic-field-driven metal-insulator phase transition in 



planar systems. Having developed a general formalism 
in the preceding sections, here we will take into account 
the effect of an external constant magnetic field jaa the 
dynamics of a spontaneous gener ation of a gap.Ea The 
general observation of Refs. 



states that in the pres- 



15 



ence of an external magnetic field, there is the generation 
of a dynamical gap connected with electron-hole pairing 
even for an arbitrary weak attraction between electrons 
and holes. Therefore in this case the gap will appear 
even if the bare coupling constant g, introduced in the 
previous section, is subcritical (in the case of a super- 
critical g, the magnetic field would enhance the already 
existent gap). This phenomenon is known as "magnetic 
catalysis" . The origin of this effect is connected with the 
dynamics of the LLL: its dynamics is effectively + 1 
(1 + 1) dimensional in 2 + 1 (3 + 1) dimensions and this 
makes the electron-hole pairing inevitable. 

Actually, this formulation is correct only in the case of 
zero temperature and zero charge density. In the pres- 
ence of temperature T and/or charge density n, there 
is a critical value of the magnetic field, B c (T,n), defin- 
ing a threshold for this effect: B has to be larger than 
B c (T,n)M While the dependence of B c {T,n) on T is 
model dependent, the value of B c (0, n) is universal for all 
values of g < g c : it is given by \e\B c = 2-Kcn/Nf anx. 
responds to the filling of the lowest Landau leveLEi 
physics of this result is quite clear: when the LLL is filled 
up, the LLL electrons are blocked and excluded from the 
pairing dynamics. In other words, in this case we loose 
the catalyst and, therefore, the effect itself. As we will 
see in Sec. VI, this point can be crucial for explaining 
the presence of an offset field B c observed in the precept 
experiments in highly oriented pyrolytic graphite.E3cJ 



A. Dynamical gap at T = and n — 

In a constant external magnetic field B, only the free 
part of the quasiparticle action is modified. In particu- 
lar, the spatial derivatives in Eq. (^|) are replaced by the 
corresponding covariant derivatives: 



cor- 
The 



d x 



c 



On 



c 



- 4 (r), 



(38) 
(39) 



where A% xt {r) = -By/2 and A e y xt (r) = Bx/2. In this 
case, the proBa.sa|tor of quasiparticlcs takes the following 
general formjl3'E3 

Git - t', r, f) = exp \-i-r ■ A ext {r') \ Git - t', f- r"), 

(40) 

Notice that while we used the symbol S for the fermion 
propagator in the case without magnetic field, we use 
the symbol G for the fermion propagator in a magnetic 
field. Let us begin by considering the propagator of free 
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quasiparticles in a magnetic field, Go(t — t', f, f). For our 
purposes, it will be useful to introduce the bare gap A& 
for these free quasiparticles. The translation invariant 
part of such a propagator, Go(t — t', f— f"), readsO 



dt d 2 k 



Gq(uj, k) = 2i exp 



2tt (2tt) 2 
f clfcl 2 



exp(— hot + ik ■ r)Go(w, k), (41) 



\eB\ J 

~ (_l)n Uurf + *b)fi{k) + 



E 



^~ A 2 -2nv 2 F \eB\/c 
In this last equation, we used the shorthand notation: 



(42) 



h{k)=V-L 




P+L n -i 




h(k) = 2v F kjL 1 n _ 1 



with the following spin projection operators: 
V± 



(43) 
(44) 

(45) 



Also, L"(z) are the generalized Laguerre polynomials. 
By definition, L n (z) = L°(z) and Lt^z) = 0. 

Let us now turn to the interactions in the presence of 
the magnetic field. In this case ike polarization effects 
could also be taken into account J23 and the modified in- 
teraction is: 



U(t,r) = S(t)- 



d 2 k exp(ifc • r) 
£o J 2tt |fc|(i + a \k\) 



e 2 nS{t) 
2e$a 



where 



a = 2iiVi 



£0^F V 



(46) 



(47) 



and the constant vq is given by 



1 



dz ( coth(z) 



3C(-0-5) 
V2vr 



1 



sinh (z) 



0.14. 



(48) 



Regarding the new notation, £(z) is the Riemann zeta 
function, H a (z) is the Struve function, and Nq(z) is 
the Bessel function of the second kind. It is noticeable 
that the instantaneous approximation for the polariza- 
tion function is justified in this case: the frequency de- 
pendence is suppressed by factors of order w/y/vF\eB\ 
(which are small in the case of the LLL dominance (see 



below)). This can be shown directly from the expression 
for the polarization function in Ref. |3^. 

Now, the gap equation for the quasiparticle propagator 
reads 



G{t,r) = G (t,r) - i I dt'd 2 ^ 



dt"d 2 ?' 

x exp — ir ■ A{r l ) — ir' ■ A(r" 

x G (i-i',f-r')7 G(<' -t",? - 
x j°G(t",^)U(^ -f")S(t' -t"). 



r") 



(49) 



The structure of this equation is essentially the same as 
in the relativistic model of Refs. |l^,^. Here, however, 
we neglect the retardation effects in the interaction po- 
tential. 

As was pointed out in Ref. [l4[ in the case of a sub- 
critical coupling constant g < g c , one should distinguish 
two different dynamical regimes. The first regime cor- 
responds to the situation with a weak coupling g, when 
it is outside the scaling region near the critical value g c . 
In this case the LLL dominates and the value of the dy- 
namical gap Ao is much less than the gap y / 2v 2 7 ,\eB\ / c 
between the Landau levels. The latter guarantees that 
the higher Landau levels decouple from the pairing dy- 
namics and the LLL dominates indeed. 

The second, strong coupling, regime is that with a 
near-critical, although subcritical, value of g. In that 
case, all Landau levels are relevant for the pairing dy- 
namics and the value of the dyn amical gap Ao is of order 
of the Landau gap y/2vp\eB\/c. 

Let us begin by considering the first regime. Then, 
the low energy dynamics is dominated by the LLL, and 
the quasiparticle propagator could be approximated as 
follows: 



G(t,f) 



i\eB\ 

47TC 



exp 



M 2 \eB\ 
4c 



il 1 ! 2 ) 



(50) 



where g(t) is unknown matrix- valued function which 
should be determined by solving the Schwinger-Dyson 
(gap) equation. By substituting the ansatz (|5^) into 
Eq. ([l9]), we derive the following equation for the Fourier 
transform of g(t): 

«r» = 9o -V) - ie 2 j fjV*(« - «V 

c\k\ 2 



d z k 



exp 



2\eB\ 



U(k). 



(51) 



The value of the bare gap is now zero in the free propaga- 
tor go(Lu). And the general structure of the function g(uj) 
is suggested by the first (LLL) term in the bare propaga- 
tor ([l2|), where now the bare gap A b should be replaced 
by the dynamical gap function A w and the wave function 
renormalization should be introduced. Thus, we have 



A^uj + A u 



A 2 ,t 



A 2 



(52) 
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One could see that the integral on the right hand side of 
Eq. (|5|) is independent of uj. This implies that A u = 1 
and the gap A w is independent of uj. By taking this into 
account, we straightforwardly derive the solution: 



_g_ v F \eB\ f 00 dkeM~k 2 ) 
° v^V c J 1 + kxo ' 



where 



Xo 



= 2V2irv gN f . 



(53) 



(54) 



In two limiting cases, xo < 1 an d Xo ^ 1, we get the 
following asymptotes: 



Ao 



g^Jv 2 F \eB\ 



Xo , Xo 



2y/2 V c \~ ^ ' 2 
(for weak coupling and small Nf) and 



A q_ h F \ eB \ m xo 

° ~ V2~V c xo 4nv N f 



eB 



(55) 



lnxo, (56) 



(for large Nf). In accordance with the general conclusion 
of Refs. HfHl in a magnetic field the gap is generated 
for any nonzero coupling constant g = e 2 / 'eqVf- 

One can sec that for a sufficiently small g = e 2 /sqVf 
in expression ( p5| ) and for a sufficiently large Nf in (56), 
the LLL approximation is sclfconsistent indeed: in both 
cases, the gap Ao can be made much less than the Lan- 
dau gap (scale) L = y / 2v F \eB\/c. We emphasize that the 
second solution (56), obtained also in Ref. H, corresponds 



to the regime with a large Nf and not to the strong cou- 
pling regime with a large g and Nf of order one. Indeed, 
taking g to be large enough in expression j5q), one gets 
the gap A exceeding the Landau scale L, i.e., for large 
g the self-consistency of the LLL dominance approxima- 
tion is lost. We will discuss the strong coupling regime 
below. 

What is the energy scale the coupling constant g re- 
lates to in this problem? It is the Landau scale L. The 
argument supporting this goes as follows. There are two, 
dynamically very different, scale regions in this problem. 
One is the region with the energy scale above the Landau 
scale L and below the ultraviolet cutoff A, defined by the 
lattice size. In that region, the dynamics is essentially the 
same as in the theory without magnetic field. In particu- 
lar, the running coupling decreases logarithmically with 
the energy scale thereJH Another is the region below the 
Landau scale. In that region, the magnetic field dramat- 
ically changes the dynamics, in particular, the behavior 
of the running coupling constant. As the analysis of this 
section shows, because of the magnetic field, the pairing 
dynamics (in the particle-hole channel) is dominated by 
the infrared region where uj < Ao- Therefore, the scale 
region above the Landau scale L completely decouples 
from the pairing dynamics in this case. This manifests 
itself in expression (|53l) for the gap: the only relevant 



scale is the Landau scale L there. Since the effect of 
the running of the coupling is taken into account by the 
polarization function in the gap equation, we conclude 
that the coupling g indeed relates to the Landau scale 
in this problem. Notice that it can be somewhat smaller 
than the bare coupling constant git) related to the scale 
t. Taking t = 2.4 eV in graphite (the width of its energy 
band) and using the equation for the running coupling 
from Ref. 7, one obtains that it is smaller by the factors 
1.2 and 1.4 than g(t) for the values of the magnetic field 
B = 10 T and B = 0.1 T, respectively. 

Now let us turn to the second, strong coupling, dynam- 
ical regime. In reduced QED, the gap equation in this 
regime includes the contributions of all the Landau levels 
and becomes very formidable. Still, one can estimate the 
value of the gap: since there are no small parameters in 
this regime for moderate values of Nf, the gap should be 
of the order of the Landau scale L. This conclusion is sup- 
ported by studying this regime in a simplcpjnodel, (2+1)- 
dimensional Nambu-Jona-Lasinio modehca In that case, 
in the critical regime, the gap equals Aq ~ 0.32L, where 
the Landau scale in that relativistic model, with v F = c, 

VWb\ 



VI, the criti- 



1S L = y/2c\eB\. As we will see in Sec. 
cal dynamical regime can be relevant for the magnetic- 
field-driven .phase transition in highly oriented pyrolytic 
graphite. EiEI 



B. Dynamical gap at T / or / 

By making use of the Matsubara formalism, it is easy 
to generalize the gap equation for the case of a finite 
temperature and a nonzero chemical potential. Without 
going into all details, let us write down the final equation, 



At(aO 



eB 



dk exp(— k 2 
1 + k\o 



sinh^ 



cosh^M 



COsh U; 



(57) 



In the LLL dominance approximation, the expression for 
charge density of carriers in terms of the chemical poten- 
tial is 



N f \eB\ 



sinh ^ 



2nc cosh 



cosh ^ 



(58) 



We assume that, in the model at hand, the charge density 
of carries (i.e., n — n e i — rih) is a fixed constant. Then, 
the expression for the chemical potential reads 



sinh - = - 



vb 



cosh 



At(m) 
T 



, sinh' 



T 



(59) 
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cosh-^ = 




where 



v B = 



2ncn 
Nf\eB\ 



Be 
B 



(60) 



(61) 



is the filling factor. 

By making use of the expression for the chemical po- 
tential in terms of the filling factor vb 1 we rewrite the 
gap equation in the following convenient form: 



A T (v B ) 



2v / 2e 



\eB\ 



dk exp(— k 2 ) 



(1 — v B ) sinh 



1 + k\o 



cosh 



T 



1 



, sinh 



2 A T Ob) 
T 



(62) 



Let us first consider the case of zero temperature. Then 
the gap equation takes a very simple form: 



A (i/fl) 



1 - v B 



e 2 (l 



A 
vb) 



2V2i 



£o 



eB 



dk exp(-/s 2 ) 
1 + k XQ 



(63) 



where Ao = Aq(0) is the value of the gap at the zero 
filling factor. Since we choose the vacuum in which the 



gap is positive (see Sec. Ill A), this equation implies that 
for Vb > 1, there is no solution with a nonzero gap, 
i.e. the symmetry is restored. The condition vb = 1 
determines the critical density, n c = Nf\eB\/2irc. The 
density n c corresponds to the filling of the LLL and, as 
was discussed at the beginning of this section, this value 
is universal for all subcritical values of the coupling con- 
stant g. Reversing the roles of n and B, one can say 
that, for a fixed value of the density n, the critical value 
of the magnetic field is \eB c \ = 2itcn/Nf: a dynamical 
gap occurs only for magnetic fields B larger than B c . 

The critical temperature is determined from Eq. j6^ ) 
with At(^b) = 0: 



4) 



4-y/2, 



eB 



dk exp(— fc 2 
1 + fcXo 



(64) 



At a fixed density n, this equation implies that, as it 
should be, the critical temperature is zero for magnetic 
fields weaker than the critical value B c determined above. 
For magnetic fields B stronger than B c , T c grows with B 
(see Fig. ^ in Sec. VI). As we will see in Sec. VI, these 
results can be important for explaining-exDerimental data 
in highly oriented pyrolytic graphite.Ea Eil 

Though here we considered only the dynamical regime 
with the LLL dominance, it is reasonable to assume that 



the qualitative picture will remain the same also in the 
case of the scaling dynamical regime, with the near- 
critical coupling constant g. This is in particular sup- 
ported by the fact of the universality of the critical value 
B c . 

Before concluding this section, let us mention that the 
the gap equation could also be rewritten in the following 
form: 



A T (u B ) 



2T C sinh 



cosh 



1 



, sinh 



2 A T (v B ) 



(65) 



where the relation ( p4| ) was taken into account. The last 
form of the gap equation will be the most convenient 
for using in numerical calculations of conductivity; see 
Sec.lVB. 



V. CONDUCTIVITY AND RESISTIVITY 

Conductivity and resistivity are major players in expe- 
rimental detecting the magnetic field driven semimetal- 
insulator phase transition in graphite JljEd In this sec- 
tion, we will calculate them in reduced QED3 + i using 
the results obtained in the previous sections. We will 
consider both cases of zero and nonzero external mag- 
netic fields. While the former case is interesting in itself, 
it will also serve us as an important reference point for 
the latter. Tne main conclusion of this section is that 
there is a clear signature of the phase transition seen in 
the behavior of the conductivity (resistivity) as a func- 
tion of temperature. More precisely, we find that 

1 . if the phase transition is of the first order, there is 
a discontinuity in the conductivity (resistivity) at 
a critical temperature T c ; 

2. if the phase transition is of the second order and 
the scaling properties are correctly described by the 
mean-field approximation, the conductivity (resis- 
tivity) exhibits a kink behavior at the critical tem- 
perature; 

3. at last, if the phase transition is a non- mean- field 
second order one, the conductivity (resistivity) is 
a smooth function at the critical temperature T c , 
while a singularity occurs in its higher derivatives 
at T = T c . 

Besides, our calculations show that in this particular 
model, the flavor phase transition, restoring the fla- 
vor symmetry U(2Np), does not look as a semimetal- 
insulator phase transition i/there is no external magnetic 
field. On the other hand, in the presence of a magnetic 
field, in many cases it does look as a semimetal-insulator 
phase transition. 
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A. No Magnetic Field 

In calculation of transport coefficients, it is very useful 
to utilize the spectral representation of the quasinarticle 
Green function. The latter is defined as follows:Eil 



S{iu) n , k) = i 



dujA(uj, k) 
iuj n — jl — LU 



(66) 



In the reduced planar model described in Sec. ||, we de- 

j°E-kj + A 7°£ + fc7-A 



A(u,k) 



2nE 



(uj-E) 2 +r 2 (uj + E) 2 +T 2 



(67) 



where E = y v 2 F k 2 + A 2 and, throughout this section, 
we use the symbol A for the gap, i.e., A = A T (fx) 
and A = At{vb) for the cases with no and with mag- 
netic field, respectively. Notice that we introduced a 
phenomenological width parameter F without which the 
calculation of conductivities would be meaningless. A 
finite width parameter appears as a result of interac- 
tions, and scattering on impurities, in particular. In 
general, the width T is defined through the fermion self- 
energy as r(a>) = —lmS R (uj). Thus, it is a frequency 
(as well as temperature and magnetic field) dependent 
quantity. Like the dynamical gap itself, it should be self- 
consistently determined from the Schwinger-Dyson equa- 
tions. At low temperatures, usually it could be modeled 
by a constant phenomenological parameter. Therefore, 
instead of considering an additional Schwinger-Dyson 
type equation, we choose a constant parameter T and 
view this as yet another approximation. 

In terms of the spectral function, the charge density of 
carriers reads 



1 f d 2 k 



n = — 
2 



(2n) 2 



duj ^tanh 



uj + (J, 
2T 



- 1 tr 



l°A(oj,k) . 
(68) 



The conductivity tensor is defined as follows: 



a tj = lim 



(69) 



where II^(f2) is the retarded current-current correlation 
function which is also given in terms of the spectral func- 
tion, 



Uij(fl + ie) 



2„,2 



e v 



cLodw' 



tanh ^# - tanh ^ 



2T 



2T 



d 2 k 



rtr 



lj — uj 1 + Q + ie 
lt A(uo,k) lo A(uo'M. (70) 



The vertex corrections were neglected in this expression. 
Formally, they are suppressed by a power of 1/Nf. Of 
course, in the case of graphite (with Nf = 2), the vertex 
contributions may nevertheless play an important role.Ej 
This question should be studied in more detail, but it 
is outside the scope of the present paper. In absence 
of a magnetic field, the conductivity tensor has only the 
diagonal components, a = a xx = a yy . Both components 
are equal as a result of rotational invariance of the model. 
The explicit expression of the conductivity, in this case, 
reads 



e 2 N f 
4ir 2 T 



T 2 duj 



dx 



{x + u 2 + T 2 ) 2 -Au 2 A 2 



C0S h 2 SH* J — [(a + lo 2 + T 2 ) 2 - Axu 2 } 2 

Z1 A 2 



8ir 2 T J cosh 

— oo 



1 



I — — arctan ■ 



2T 

r 2 + A 2 - 



2 \u> T 



2Ur 



(71) 



where T is the width parameter, and the density of car- 
riers is defined by the following relation: 



2ir 2 v 2 F 



du 



T 2 



dEE 



A 



tanh^ + ^+tanh' + f ^ 



2T 



2T 



(72) 



In the limit r — * 0, these two expressions reduce to 

- A 2 

-^6(cu 2 - A 2 ) 



e 2 N f f dw 



167T7T 



cosh 2 age 



e 2 N f f dco 



16ttTF J oj 

A 



cosh 2 ^ 



(73) 



and 



2irv 2 F 



dEE 



tanh 



N f T 2 sinh f 



fi + E 
2T 

dxx 



tanh 



H- E 
2T 



cosh x + cosh if 



nvi 



A l + cxp(^) 
T n i + exp(-^) 



Li 



Li 2 ( -e*T 



(74) 



where lA-z^z) is the dilogarithm function. As one could 
see from the above formulas, the conductivity grows lin- 
early with temperature when the temperature is large, 
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e 2 N f T 

4?r r 



In 2, for T 



(75) 



Notice, however, that the expression for the conductiv- 
ity in Eq. d73]), derived for the r — > case, fails when 
temperatures become very small. The correct result for 
small temperatures could be derived from Eq. (|n]). It 
reads 



2tt 2 



1 



A 2 + r 2 



x | arctan 

2 



2| M |r 
r 2 + a 2 - f 



2|/i|r 

The density in that same limit is 
- N f 



2nvj, 



- A>gn(^(^ 



A 



(76) 



(77) 



The interplay between the density of carriers and the 
width r is characterized by the following dimensionless 
parameter: 



(78) 



In the two opposite limits of a clean or dirty system, this 
parameter is either large or small, respectively. Then, 
the corresponding zero temperature asymptotes for the 
conductivity take the form: 




e 2 N f r] I 
4tt y 2nv 2 F n + N f A 2 ' 



2irVpn 



for r) ^> 1, and 
e 2 N f 



r /tt r 

1 H — arctan — — 

2A \ 2 2A 



(79) 



(80) 



for r\ <C 1. The last expression was derived under the 
assumption that n <C A 2 /v F . Finally, in the strict limit 
of zero density (i.e., /i = 0), we derive 
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It should be emphasized that the strict case of zero den- 
sity corresponds to [i — (rather than \i = A as might 
be suggested by taking the limit T — > first, and then 
n — * 0). To understand this better, one should look at 
the temperature dependence of the chemical potential at 
a given fixed value of the density. In particular, when 
the density of carriers is very small, the chemical poten- 
tial as a function of temperature sharply falls from its 
value /i = A at T — almost down to zero in a very 
small region of temperatures. Afterwards, it starts to 
grow. When the density gets vanishingly small, the be- 
fore mentioned region of temperatures where the chem- 
ical potential drops shrinks to zero. Thus, by making 



use of continuity argument, it is clear that the value of 
the chemical potential is zero in the limit T — > if the 
density of carriers is zero. 

It is noticeable that our r csult in Eq. @ is m agree- 
ment with the Wiedemann- Franz law, i.e., 
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(82) 



where we use the value of the thermal conductivity k 
calculated in Ref. [|| 

The numerical results for the temperature dependence 
of the conductivity are shown in Fig. ^| in the case of 
the zero gap (bold solid line), and nonzero gaps (other 
lines correspond to different values of T c ). Notice that the 
model at hand reveals an "insulator" (i.e., increasing with 
temperature) type behavior of conductivity even in the 
case of a finite density of carriers. This type of behavior 
is the consequence of using a constant value of the width 
parameter T in our model. Then, the growth of con- 
ductivity with increasing temperature is directly related 
to the increasing number of thermally excited quasipar- 
ticles. In realistic systems, of course, the width (which 
is related to the inverse scattering time) would normally 
start to grow with temperature too. In general, one might 
choose the width as function of energy and temperature, 



r(^,T) = r 



i 



T(u,Ty 



(83) 



where Tq is the zero temperature width due to impurities, 
and the other term is due to the thermal contribution. 
In this paper, for the sake of simplicity, we consider the 
simplest model with a fixed constant value of the width 
parameter. The analysis, however, could be easily gener- 
alized for any phcnomenologically motivated dependen- 
cies like that in Eq. (p^). 

In order to calculate the conductivities in the case of 
nonzero dynamical gaps, we used the gap equation ( |65| ) 
in which the critical temperature T c was treated as a free 
parameter. 

The results for the temperature dependence of the re- 
sistivity are plotted in Fig. |J. 

As one can see in Fig. || and Fig. ||, there is a kink at the 
critical point T = T c in the conductivity (resistivity) . Its 
occurrence is directly related to the mean field behavior 
of the gap A in the vicinity of the critical point, i.e., 
A - y/T c - T. Indeed, as follows from Eq. (JtTJ) , the 
conductivity a depends on A 2 and there is a linear in 
A 2 term in it as A 2 — > 0. Therefore, its derivative with 
respect to temperature has a finite discontinuity at the 
critical point T = T c . 
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FIG. 3. The conductivity as a function of temperature for 
the zero magnetic field case. The bold solid line corresponds 
to the case without a gap. The other lines correspond to 
nonzero dynamical gaps and different values of T c . Conduc- 
tivity is measured in units of e 2 , both temperature and width 
r are measured in Kelvin, and density n is measured in cm -2 . 
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FIG. 4. The resistivity as a function of temperature for 
the zero magnetic field case. The bold solid line corresponds 
to the case without a gap. The other lines correspond to 
nonzero dynamical gaps and different values of T c . Resistivity 
is measured in units of e -2 , both temperature and width Y 
are measured in Kelvin, and density n is measured in cm -2 . 



The mean field behavior may change if higher order, 
1/Nf, corrections (fluctuations) are taken into account. 
The fluctuations could either change the phase transi- 
tion to a first order one, with a discontinuity in A at the 
phase transition point, or to a non-mean-field second or- 
der phase transition, with the scaling law A ~ (T c — T) v 
where v > 1/2. While in the former case a discontinuity 
will appear in the conductivity (resistivity), in the latter 
case the conductivity (resistivity) will be a smooth func- 
tion of temperature, and a singularity will move to its 
higher derivatives. 



Another noticeable point is that in the case with no 
magnetic field, the flavor phase transition does not look 
as a semimetal-insulator one. Indeed, as one can see in 
in Fig. H and Fig. ||, an "insulator" type behavior does 
not change at the critical point. As we will see in the 
next subsections, the occurrence of a magnetic field will 
drastically change this feature of the phase transition. 

Before concluding this section, let us also mention that 
the conductivity (as well as the resistivity) become more 
sensitive to the appearance of a dynamical gap when the 
density of carriers decreases. To support this statement, 
we plotted the conductivity and resistivity for two dif- 
ferent finite values of carrier densities which differ by a 
factor 10, see lower parts of Figs. || and || 
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As we see from Fig. |], the temperature dependence 
of the resistivity develops a minimum when the value of 
the gap is sufficiently large. Comparing this temperature 
dependence with the experimental data, we might even 
suggest that the studied graphite samples are better de- 
scribed by the model with a nonzero dynamical gap even 
in absence of a magnetic field. The effect of an external 
field is studied in the following subsections. 



B. Conductivity tensor. Nonzero magnetic field. 

Let us now turn to the analysis of the conductivity in 
the case with an external magnetic field. The spectral 
function A(u>, k) of the translation invariant part of the 
quasiparticle propagator in a magnetic field [see Eqs. flU] ) 
and (|42|)1 is given by 
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(w - m„) 2 + r 2 

( 7 °M n -A)/ 1 (fc)-/ 2 (fc) 
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(84) 



where M n = -y/A 2 + 2nv' F \eB\/c and the functions fi(k) 

and f2(k) were defined earlier in Eqs. ( f43| ) and (Q). In 
an external magnetic field, the conductivity is a tensor 
quantity. The diagonal and off-diagonal components of 
conductivity read 
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respectively. Here the parameter T gives the energy 
width of Landau levels, and the filling factor vb [re- 
lated to the density of carriers by the relation: vb = 
2ncn/(Nf\eB\)] is defined as follows: 



VB 



— oo 
oo 



duo a + lj 
— tanh — — — 
2tt 2T 



{ V(^-M„) 2 +L 2 



(w - A) 2 + L 2 
+ (w — * — w) 



+ (id — * — Ld) 

(87) 



The sum over Landau levels in Eq. (j8q ) could be per- 
formed explicitly, and the result is given in terms of the 
digamma function, ip(z), as follows: 
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The high temperature asy mptote that follows from the 
representation in Eq. ( p8[ ) is the same as in the case of 
zero magnetic field, given in Eq. ( |75| ) . The limit T — > is 
different from that in Eq. (|76|). It is given by the following 
expression: 
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(89) 



It is interesting to note, however, that for zero value of 
the gap and zero density of carriers (i.e., A = and 
|U = 0), this last expression becomes identical with the 
expression for the conductivity in absence of a magnetic 
field given in Eq. Q8l|) . 

In the limit of narrow width, F — > 0, the above expres- 
sions reduce down to 
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(91) 



for the filling factor. 

In order to understand the effect of a dynamical gap 
on the behavior of conductivity as a function of temper- 
ature, it is helpful to start from the case of a vanishing 
density of carriers (i.e., vb — 0). When vb = 0, the Hall 
conductivity is absent, and the resistivity is determined 
by (J xx component alone. The plot of the conductivity as 
a function of temperature is given in Fig. [| 
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FIG. 5. The diagonal component of conductivity as a func- 
tion of temperature for zero value of carrier density and a 
nonzero magnetic field, B — 8 Tesla. The bold solid line cor- 
responds to the case without a gap. Other lines correspond 
to nonzero dynamical gaps and different values of T c . Con- 
ductivity is measured in units of e 2 , both temperature and 
width F are measured in Kelvin. 

The bold solid line corresponds to the case without a 
dynamical gap, while the other lines correspond to dif- 
ferent choices of the gap magnitude. In the numerical 
analysis, we used the gap equation in Eq. (|65|), keeping 
the value of T c as a free parameter. 

When there are no free carriers, the low temperature 
dependence of the diagonal component of conductivity 
is very sensitive to the presence of a gap. In absence 
of a gap, the conductivity becomes infinitely large when 
T — > 0. At the same time, it is zero in the same limit 
when there is even an arbitrarily small gap A. 

An important fact is that, unlike the case without 
magnetic field, the conductivity exhibits a (semi)mctallic 
type behavior for zero gap and not too high temperatures 
(T < 200 K in Fig. ||). As a result, for not too large values 
of the critical temperature T c (T c < 200 K in Fig. |), the 
flavor phase transition looks as a conventional semimetal- 
insulator one, when the insulator type behavior below T c 
(nonzero gap) is replaced by the metallic type in a range 
of temperatures just above T c (zero gap) (see Fig. ||). 

A typical conductivity for a nonzero value of the fill- 
ing factor ub (i.e., nonzero charge density) is shown in 
Fig. |6[ In this case, the behaviors of the conductivity 
for a nonzero gap and zero gap are more similar than in 
the case of vb — (compare with Fig. ||). The presence 
of a gap, however, can substantially reduce the value of 
conductivity in the whole range of temperatures below 
T c . It is important that, like in the case with vb = 0, 
the flavor phase transition looks as a semimetal-insulator 
one for not too large values of the critical temperature T c 
(see Fig. ^j. It is the most important conclusion of this 
subsection. 



The same arguments as in the end of Sec. VA show 
that the occurrence of the kink in the conductivity at 
T — T c reflects the mean-field behavior of the gap in 
the vicinity of the critical point, A ~ \/T c — T. The 
1 /Nf fluctuations may change the character of the phase 
transition, leading either to a discontinuity in the con- 



ductivity a(T) at T = T c (a first order phase transition) 
or to a smooth function u(T) (a non-mean-field second 
order phase transition). 
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FIG. 6. The diagonal component of conductivity as a func- 
tion of temperature for two different densities and a nonzero 
magnetic field, B — 8 Tesla. The bold solid line corresponds 
to the case without a gap. Other lines correspond to nonzero 
dynamical gaps and different values of T c . Conductivity is 
measured in units of e 2 , both temperature and width V are 
measured in Kelvin, and density n is measured in cm -2 . 



C. Resistivity tensor 

In this subsection, we study the temperature depen- 
dence of the resistivity. 

In terms of conductivities, the diagonal component of 
the resistivity reads 



Px 



(92) 



In order to understand the general behavior of the resis- 
tivity, below we perform a set of numerical calculations. 

Before presenting the results, it is instructive to notice 
that there exist two opposite regimes of dynamics con- 
trolled by the value of the charge density. In particular, 
at small density, when the Hall conductivity a xy is negli- 
gible compared to a xx , the resistivity in Eq. (^2| behaves 
as l/a xx . On the other hand, at sufficiently large density, 
when the Hall conductivity dominates over the diagonal 



component, the resistivity p xx « o xxl 
that the Hall conductivity [see Eq. 



By recalling 
])] is independent 
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of temperature, the general features of the temperature 
dependence of p xx will be the same as of 1/ a xx and a xx 
in the mentioned two regimes, respectively. 

Now, let us present the numerical results. We begin by 
considering the case of zero density. In this case the Hall 
conductivity equals zero and the resistivity p xx equals 
1/cTxx- The temperature dependence of the resistivity is 
shown in Fig. |^ (compare with Fig. |^). 
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FIG. 7. The resistivity as a function of temperature for zero 
value of carrier density and a nonzero magnetic field, 5 = 8 
Tesla. The bold solid line corresponds to the case without a 
gap. Other lines correspond to nonzero dynamical gaps and 
different values of T c . Resistivity is measured in units of e -2 , 
both temperature and width V are measured in Kelvin. 

The bold solid line corresponds to a model with the 
vanishing dynamical gap in the quasiparticle spectrum. 
The other curves correspond to three different choices of 
the dynamical gap. The temperature dependence of the 
dynamical gap is given by the gap equation ( |65| ) where 
the value of T c was treated as a phenomenological pa- 
rameter. 

As one can see in Fig. |7| for zero dynamical gap (the 
bold solid line) the resistivity has a metallic type behavior 
for not too high temperatures (T < 0.2vF\/\e-B\/c) and 
an insulator type behavior at high temperatures (T > 
0.2vF\/\eB\/c). This type of temperature dependence 
is driven by the magnetic field alone and is not related 
to the generation of a dynamical gap. Such a crossover 
from the metallic type behavior (low temperatures) to 
the insulator one (high temperatures) can be easily dis- 
tinguished from the flavor phase transition taking place 
at not too high T c . Indeed, when T c < 0.2vpy/ \eB\/c, 
one can see from Fig. [t] that it corresponds to the op- 
posite, conventional, transition, when the insulator type 
behavior below T c (nonzero gap) is replaced by the metal- 
lic type in a range of temperatures just above T c (zero 
gap). 

Let us now proceed to the case of a nonzero charge 
density. As we mentioned at the beginning of this sub- 
section, there are two different regimes that appear in 
the limits of small and large densities, respectively. Our 
results in Fig. || illustrate these regimes. 
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FIG. 8. The diagonal component of resistivity as a function 
of temperature for four different densities and a nonzero mag- 
netic field, B = 8 Tesla. The bold solid line corresponds to the 
case without a gap. Other lines correspond to nonzero dynam- 
ical gaps and different values of T c . Resistivity is measured 
in units of e -2 , both temperature and width F are measured 
in Kelvin, and density n is measured in cm -2 . 

At small density (see the upper panel in Fig. ||), the 
resistivity behaves as l/a xx almost at all temperatures 
(compare with Fig. |^), except for a finite region where 
a xx becomes very small due to the generation of a dy- 
namical gap. In this region, even a small value of the Hall 
conductivity could dominate over a xx . This is seen as the 
appearance of a local minimum (between two maxima) in 
the temperature dependence of the resistivity. Of course, 
in the absence of a dynamical gap (see the solid line in the 
upper panel in Fig. |^), the resistivity remains essentially 
the same as that at n — (compare with Fig. |7|). 

This picture changes dramatically with increasing the 
density. As one can see in the two lower panels in Fig. 
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with increasing the density the resistivity gradually ap- 
proaches the regime where p xx ~ o~ xx . 

For critical temperatures T c < 0.2vjry/]eB]/c (i.e., for 
dependencies represented by thin solid and dashed lines) , 
these two regimes correspond to two essentially different 
metal-insulator phase transitions. In the case of a small 
density (see the upper panel in Fig. ||), it is a conven- 
tional phase transition with insulator type and metallic 
type behaviors at temperatures below and just above T c , 
respectively. On the other hand, for a large density (see 
the lower panel in Fig. |8|), the "inverse" metal-insulator 
phase transition (with metal type dependence just below 
T c and insulator type just above T c ) is realized. 

Therefore we conclude that, in the presence of a mag- 
netic field, a dynamical gap in the quasiparticle spec- 
trum can indeed lead to a change of the insulator type 
dependence of p xx (T) to the metallic one. However, the 
nonzero Hall conductivity at finite density n complicates 
the picture and can lead to different types of metal- 
insulator phase transitions for small and large values of 
the charge density. 

When the metal-insulator phase transition is truly a 
mean-field one, its clear signature is a kink in the resis- 
tivity p X x(T) at the critical point T = T c . As has been 
already pointed out in the previous subsections, the 1/Nf 
fluctuations can change this feature, leading either to a 
discontinuity in the resistivity at T — T c (a first order 
phase transition) or to a smooth function p xx (T) (a non- 
mean- field continuous phase transition). 



VI. METAL-INSULATOR PHASE TRANSITION 
IN HIGHLY ORIENTED PYROLYTIC 
GRAPHITE 

The main motivation of this study was the experimen- 
tal data reported in Refs. ^9|-^l]. It was observed that 
samples of highly oriented pyrolytic graphite in an ex- 
ternal magnetic field show a qualitative change of their 
resistivity as a function of temperature, that was inter- 
preted as a metal-insulator phase transition. The effect is 
clearly seen only for a magnetic field perpendicular to the 
basal plane, suggesting that the orbital motion of quasi- 
particles is responsible for the change of the conductivity 
dependence. 

In this section we will attempt to explain qualitatively 
the main features of the above mentioned experimen- 
tal data in the light of the magnetic catalysis idea. We 
should note that the first step in this direction was made 
in Ref. [|. Here we go into further details utilizing the 
rather complete description of the magnetic catalysis in 
planar systems and its effect on the temperature depen- 
dence of their conductivity and resistivity given in the 
previous sections. 

First of all, the analysis made in Sec. |y| shows that, in 
the presence of a magnetic field, the flavor phase tran- 
sition in planar systems can indeed manifest itself as a 
metal-insulator phase transition in the behavior of their 



resistivity p(T) as a function of temperature. A notice- 
able fact is the existence of clearly distinguishable sig- 
natures of different types of the phase transition: the 
presence of a discontinuity and a kink in the resistivity 
p(T) at the critical point T = T C in the cases of first or- 
der and mean-field phase transitions, respectively, and a 
smooth behavior of p(T) at T = T c for a non-mean-field 
continuous phase transition. To the best of our knowl- 
edge, so far there have been no experiments reporting 
observations of a singular behavior of p(T) at the criti- 
cal point. At this stage, however, it would be premature 
to conclude that the observed phase transition is a con- 
tinuous non-mean-field one. This point deserves further 
experimental study. 

A very interesting experimental observation made in 
Refs. pj|-pT| (and, to the best of our knowledge, has not 
been explained) is the existence of a finite "offset" mag- 
netic field B c . The value B c determines the threshold 
B = B c for a qualitative change of the resistivity at zero 
temperature. More precisely, based on the experimen- 
tal data, it was revealedRil that the approximate rela- 
tion for the critical temperature as a function of B reads 
T C (B) ~ \JB — B c . This relation implies that at zero 
temperature the phase transition happens only when the 
magnetic field exceeds the threshold value B = B c . 

It is remarkable that, as was emphasized in Sec. IV, 
the existence of such a threshold B c is a robust conse- 
quence of the mechanism of the magnetic catalysis. As 
was pointed out there, the value B c is directly related to 
a nonzero charge density n of carriers, 



\eB c \ = 



2ncn 

w 



(93) 



and this relation is exact. For example, by taking B c — 
2.6 x 10 4 G, which was obtained in one of the experiments 
as an (upper) estimate of a critical value above which the 
generation of a gap presumably occurred, t£l we derive the 
corresponding charge density (we use Nf = 2) 



\eB r 



TTC 



1.25 x lO -5 ^ -2 . 



(94) 



This should be compared with the charge density (per 
unit area of a layer) of carriers in the used sample of 
graphite. By noting that the area per carbon atom 
in a layer is S — v / 3a 2 /4 where the lattice spacing is 

a sa 2.46A,E1 we conclude that the density in Eq. Q) 
corresponds to n ~ 3.3 x 10~ 5 units of charge per atom. 
While we do not know the exact densijtjj of the sample 
used, the given estimate is not unlikelyHfl 

Notice that the relation T C (B) ~ ^B - B c = 
"v/1 — vbVb, used in Ref. ^l], qualitatively differs from 
our Eq. (|64|). It is quite remarkable, however, that the 
dependence T C (B) in Eq. ( |64| ) is nearly the same numer- 
ically as the simple square root relation suggested by the 
experimental data, see Fig. |9[ 
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FIG. 9. The critical temperature as a function of a mag- 
netic field. The solid line gives the dependence in Eq. (p4j), 
while the dashed line corresponds to a dependence with a sim- 
ple field offset, i.e., T c ~ \fB — B c . To plot the figure we used 
the following parameters: Nf = 2, c/vf = 375, £o = 2.4 and 
B c = 2.6 x 10 4 G. 

Since the relation \eB c \ — 2ircn/Nf is exact in the 
dynamics of the magnetic catalysis, its experimental ver- 
ification would be a critical check of the validity of the 
scenario of the magnetic catalysis in highly oriented py- 
rolytic graphite. 

Another noticeable experimental observation is that 
the scale of the critical temperature is set by the en- 
ergy distance between the Landau levels (the Landau 
scale) Jiff til Therefore, if the underlying physics of the 
transition is related to a dynamical generation of a gap, 
the typical values of the gap should also be of the same 
order as the Landau scale. As we discussed at length in 
Sec. IV , for the mechanism of the the magnetic catalysis 
this fact implies that the pairing dynamics corresponds to 
the strong coupling regime. This in turn implies that all 
(or many) Landau levels determine the pairing dynamics 
in this case. In connection with that, we would like to 
point out that, as the numerical analysis done in Sec. |v| 
shows, the contribution of higher Landau levels into the 
conductivity and resistivity become indeed important for 
values of the critical temperature T c of the order of the 
Landau scale. 

There still remain some unresolved issues in the inter- 
pretation of the experimental datajniiighly oriented py- 
rolytic graphite in a magnetic field0~0 The most impor- 
tant of them is the observation of weak ferromagnetism in 
that system (for some speculations concerning its origin 
see Ref. H). We hope to consider this issue elsewhere. 



VII. CONCLUSION 



systems was described in detail. This allowed us to con- 
clude that, in the presence of a magnetic field, the flavor 
phase transition in planar systems can indeed manifest 
itself as a metal-insulator phase transition in the behav- 
ior of the resistivity p(T, B) as a function of the magnetic 
field and temperature. It was also shown that there exist 
clearly distinguishable signatures of different types of the 
phase transition. While the resistivity p(T) is a smooth 
function at the critical point T = T c in the case of a 
non-mean-field continuous phase transition, there are a 
discontinuity and a kink in p(T) at T = T c in the cases of 
the first order and mean-field phase transitions, respec- 
tively, nn 

Based on the experimental data,ta~Ell it has been re- 
cently argued that highly oriented pyrolytic graphite 
shows up a metal-insulator phase transition, driven by an 
external magnetic field.H This might be a nonrelativistic 
realization of the phenomenon of the magnetic catalysis 
originally established in Refs. 14 ll] in relativistic sys- 
tems. In this paper we studied this possibility rather 
in detail, elaborating the theory of the magnetic catal- 
ysis in nonrelativistic planar systems and analyzing the 
temperature behavior of the resistivity (conductivity) in 
these systems. The conclusion of the present analysis 
concerning the possibility of this scenario in highly ori- 
ented pyrolytic graphite is quite positive. 

One of the central results of this paper is establishing 
the exact relation (|93| ) for the critical (threshold) value of 
the magnetic field at zero temperature in these systems. 
An experimental verification of this result would be a 
crucial test for the present theory. 

Another conclusion of our investigation is that a 
nonzero magnetic field alone (even without producing a 
dynamical gap) can drastically change the general be- 
havior of the resistivity as a function of temperature. In 
particular, in our simplest model with a constant value 
of the width parameter, the semiconductor type depen- 
dence of the resistivity [i.e., p(T) decreasing with increas- 
ing temperature] , seen in the absence of a magnetic field, 
can be replaced by a metallic type behavior [i.e., p(T) 
increasing with temperature] in the region of not too 
high temperatures, when a nonzero field is turned on. 
In fact, at zero charge density of carriers, this change of 
behavior always happens in the range of temperatures 
< T < 0.2vF\J\eB\/c. This is also seen at finite but 
small densities when the diagonal component of conduc- 
tivity dominates over the Hall conductivity. 

We expect that the results of this paper will be useful 
for a wide class of condensed matter planar systems. 



In this paper we developed a theory of the magnetic- 
field-driven metal-insulator phase transition in planar 
systems, based on reduced QED. The general structure of 
the phase diagram of such systems was established in two 
cases, with and without an external magnetic field. The 
behavior of the electric conductivity (resistivity) in these 
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APPENDIX A: SYMMETRY OF 
(2 + l)-DIMENSIONAL FERMIONS 

In this Appendix we will consider the symmetry of 4- 
component fermions on a plane which carry the flavor 
index i = 1,2, ... , Nf. The three 4x4 7-matrices in 
Eq. (||) can be taken to be 



0":s 
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(Al) 


io\ 





(A2) 


102 



V 

— 1(72 J 


(A3) 



7 



Recall that in 2 + 1 dimensions, two sets of matrices 
(03, io\, 102) and (— CT3, — ioi, — 102) make inequivalent 
representations of the Clifford (Dirac) algebra 



where fx,u = 0, 1, 2 and = diag(l, -1, -1). 
There are two matrices, 



7 3 = i( ° 1 
* 1 1 



7 1 1 -1 



(A4) 



(A5) 



that anticommute with 7 , 7 1 and 7 2 . Therefore for each 
four-component spinor, there is a global U(2) symmetry 
with the generators 



I, -7 3 , 7 5 , and ^[7 3 ,7 5 ] 
% 2 



(A6) 



Since there are Nf fermion flavors, the full symmetry of 
the action (|Io| ) is U{2Nj) with the generators 

where A Q /2, with a = 0, 1, ... , N 2 — 1, are AT 2 generators 
of U(Nf). ' " 

Adding a mass (gap) term A tjjtp into the action (nO) 
would reduce the U(2Nf) symmetry down to the U (Nfjx 
U(Nf) with the generators 

A" A Q 1 



2 2 



[7 3 ,7 5 L 



(A8) 



with a = 0, 1, ... , Nf — 1. This implies that the dynam- 
ical generation of the fermion gap leads to the sponta- 
neous breakdown of the U(2Nf) down to the U(Nf) x 
U(N f ). 



APPENDIX B: DERIVATION OF 
POLARIZATION FUNCTION AND GAP 
EQUATION 

In this Appendix, we give the details of the calculations 
of the time component of the gauge field polarization 
function, as well as the derivation of the gap equation 
at finite chemical potential and finite temperature. We 
will consider only the case of zero magnetic field. The 
polarization function and the gap equation in (2 + 1)- 
dimensional QED with an external magnetic field were 
given in Ref. B(J where the method of Ref. [lq was used. 



1. Polarization function 



The general expression of the time component of the 
vacuum polarization function is given byEj 



2tt o 

U(n m ,p) = — e 2 

£0 



n=~oo J v 

y S(Q m + uj„, p + k)~f S (u> n , k) 



d 2 k 



x tr 



(Bl) 



where S(uj n ,k) is the fermionic quasiparticle propagator 
whose explicit form reads 



S(uj n ,k) 



(iuj n - m)7o + (k ■ 7) + At (a*) 



(B2) 



In Eq. (Bl), the Matsubara frequencies are denoted by 
u> n = (2n + 1)ttT and Q, m = 2mirT. Also notice that the 
expression on the right hand side is multiplied by an ad- 
ditional factor 2ir/eo, in accordance with our definition 
of the polarization function. After taking the trace over 
the Dirac indices and using the Feynman parametriza- 
tion, we obtain 



n= — oo u 



d 2 k_ 



K + in) 2 + v 2 F p 2 x(l - x) + v 2 F k 2 + A|(/i) 

2{v 2 F k 2 + A 2 T (ii)} 

[{uj n + in) 2 + v 2 F p 2 x(l -x) + v 2 F k 2 + A 2 ,^)} 2 

By calculating the sum over n, we get 
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11(0, p) 



e 2 N f 
2e 



dx 
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v F p x{\ - x) 
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x tanh 



Y + n v 2 F k 2 + A F ([i) 



2T 



2Tcosh 2 ^ 
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(B4) 



where Y = ^/v F k 2 + v F p 2 x(l — x) + A 2 ^(fi). By chang- 
ing the integration variable, k 2 Y, and integrating by 
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parts, we finally arrive at the following convenient repre- 
sentation: 



11(0, p) = 



2Te 2 N f 



£qv f Jo 



d.r 



In 2 cosh 



Rx + fi 



^■tanh^ 



2TR X 



2T 



2T 

+ -> -/■*) 



(B5) 



where R x = ^Jv 2 F p 2 x{\ — x) + A F ([i). 



2. Gap equation 

The general Schwinger-Dyson (gap) equation for the 
quasiparticle propagator reads 



T 



S 1 (uj m ,p) = S 1 (uj m ,p) 

x 7 5K,fc)7°C/(p-fc) 



n=— oo 



d 2 k 



(B6) 



By neglecting the wave function rcnormalization,il 
derive the following gap equation: 



wo 



A(p) = 



e 2 T 



+ 00 



27T£o ^ 
n —-co 



A(k)d 2 k 



(u n + i^) 2 + v 2 F k 2 + Af,(/x) 



\p-k\ + Tl(0,$-k)' 



(B7) 



where Af(/i) = A(p)\ p= o- Here the interaction is taken 
in the so-called instantaneous exchange approximation. 
This means that the retardation effects of the gauge field 
are neglected which is justified in a nonrelativistic model. 

By neglecting the dependence of the gap on the Mat- 
subara frequency, we could perform the sum over n ex- 
plicitly. Then, the result reads 



sinh 



vre^ r d 2 k A{k) 

{P) ~ e J (2ir) 2 E k cosh + cosh f 
1 



\p- k\ +U(0,p- k) 



(B8) 



where Ek — y 'v 2 F k 2 + Aj,(fi). By using the standard 
approximation for the kernel of the integral equation, 
f(\p- fc|) -> /(p)6»(p - fe) + /(fc)6»(fc - p), we obtain the 
following gap equation: 



A( P ) 



2e v F 
6{p - k) 



dkA{k) 



sinh^ 



cosh ^ + cosh f 



6{k-p) 



p + n(0,p) fc + n(0,fc) 



(B9) 



where the infrared cutoff e is given by a larger value of 
At(h)/vf or -\//i 2 — A\{h)/vf, and where we also uti- 
lized the bifurcation method in which a nonlinear gap 



equation is replaced by a linear a pprox imation (compare 
with the discussion in subsection III B ). This is achieved 
by substituting the trivial value of the gap in Ek and in- 
troducing an infr ared cutoff in the integral on the right 
hand side of Eq. @. 



APPENDIX C: DERIVATION OF EFFECTIVE 
POTENTIAL AT /i/0 

In this appendix, we will construct the effective po- 
tential of the composite field a = —(i/ji/j) by using the 
method of Ref. |j| For the purposes of this paper, it is 
sufficient to consider only the case of a nonzero chemical 
potential. The generalization to some other cases (for ex- 
ample, with an external magnetic field) is also possible, 
see for example Ref. |36|. 

In order to derive the effective potential as a function of 
the composite field er, one should introduce a term with a 
constant external source J coupled to the corresponding 
composite operator in the action, and construct the gen- 
erating functional W(J). The effective potential, thee, 
is defined through the Legendre transform as follows SB 



V((j) = -w(J) + Jcr = / daJ(a), 



(CI) 



where a = dw(J)/dJ, w(J) = W(J) /V2+1, and V2+1 is 
the space-time volume. In the last expression, the source 
J should be regarded as a function of the field a. 

The effect of the external source J could be easily taken 
into account in the gap equation (^9|) : one should simply 
replace A p — > A p — J on the left hand side of the equation. 
Then, the solution to the equation, satisfying the infrared 
boundary condition, takes the following form: 



A /e 

' - sin 



sin Syp 



2 e 



(C2) 



where v = \J A\ — 1, e = max{A/t>i?, \/ [i 2 — A 2 /vf}, 
and 5 = arctanz/. The overall normalization of the above 
solution is fixed by choosing A p=e = A. The ultraviolet 
boundary condition, 



J= (A p+P A' p )\ p=A , 
on the other hand, produces the relation: 



sin(2(5) V A 



sm 



2 e 



2S 
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As it should be, the equation for the dynamical gap Ao 
is obtained in the limit of vanishing source J = 0. At 
zero chemical potential, in particular, the equation for 
the gap takes the form: 



v Avf 
x In — ; — 
2 A 



2(5. 



(C5) 
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For the derivation of the effective potential, we also need 
to know the expression for the field a. By definition, it is 
equal to the trace of the fermion propagator. Thus, we 
get 



a = -(ipi/j) = - 



N 



J 



nXvi 



Nf AVeA 
ttXvf sin(2<5) 



sm 



2 e 



(C6) 



Now, by making use of Eq. (C5), we trade the cutoff 



parameter A for Ao. After that, we derive the follow- 
ing approximate relations for the case of small v we are 
interested in: 



4 V A cvf 



a(A) 



ttvf 



4 -In 



_Ao_ 



(C7) 
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As will become clear in a moment, these two expressions 
contain all the information needed for reconstructing the 
potential. Indeed, the definition of the effective potential 
in Eq. ypw can be rewritten as follows: 



V(a) 



dA^J(A) + /(/*), 



(C9) 



where the most general integration constant f(ti) was 
added on the right hand side. This new representation 
leads to the final result, 
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+ h(/i), (cii) 



for A > \fi\/\/2. In these equations we used the fre edom 
of ch oosin g the integration constants in expressions ( C10| ) 
and (|Cll| ) as follows: 
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This choice insures that the potential is continuous at the 
matching point A = |//|/v2, and that it is normalized so 
that its partial derivative with respect to the chemical 
potential at the global minimum is equal (up to a sign) 
to the charge density: 



for < fi c , and 

dV(0,fJ.) 
<9/i 



2ttw 2 



sgn( M )0(| M | - Ao), (C14) 



S Sgn(M) ' 



(C15) 



for > fi c . Here we used the expression for the charge 
density in Eq. @. 

Now, the effective potential as a function of the c ompo - 
site fi eld a is defined parametrically through Eqs. (CIO), 
TTI) and (T 



). This dependence is graphically shown 
for a few different values of the chemical po- 
As is clear from the figure, the presence of a 



in Fig. 
tential. 

nonzero chemical potential considerably changes the be- 
havior of the effective potential. In particular, a new local 
minimum develops at the origin and its depth gradually 
increases with /i. The competition of the two minima, 
located at a = and o§ = <j(Ao), results in a first order 
phase transition. Such a transition happens when the 
depths of effective potential at its two minima become 
equal. By making use of this criterion, we derive the an- 
alytical expression for the critical value of the chemical 
potential, 



An 



(2 - V2)!/3 



1.195A . 
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